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O ■ Abstract 
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By employing monomer-resolved computer simulations and analytical considerations based 
^ I on polymer scaling theory, we analyze the conformations and interactions of multiarm star 

_^_' polymers strongly adsorbed on a smooth, two-dimensional plane. We find a stronger stretch- 

ing of the arms as well as a stronger repulsive, effective interaction than in the three dimen- 
sional case. In particular, the stai^ size scales with the number of ai^ms / as ~ /''^ and the 
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effective interaction as ~ /^, as opposed to ~ f^l^ and ~ /^/^, respectively, in three dimen- 
sions. Our results demonstrate the dramatic effect that geometric confinement can have on the 
effective interactions and the subsequent correlations of soft colloids in general, for which the 
conformation can be altered as a result of geometrical constraints imposed on them. 



Introduction 

Physical confinement of complex fluids in small pores, films, or tubes is known to induce sig- 
nificant changes in their structural and dynamical properties and phase behavior, including the 
emergence of surface-related phase transitions-^*^ the modification of bulk critical phenomena,-"^ 
novel types of ordering-"— as well as a shifting of dynamical arrest lines.— As such, the de- 
tailed microscopic understanding of this confinement-based phenomenology is of great importance 
both scientifically and technologically. Already when hard, nondeformable colloids are restricted 
within a slit'^ or on the interface between two immiscible fluids,-^ important new mechanisms 
come forward which affect the effective interactions between the constituent particles. Common 
examples are the softening of the interactions emerging from averaging-out the transversal degrees 
of freedom-ii^ or the development of strong electric dipoles due to dielectric mismatch between 
the immiscible fluids forming the interface.-!^ To date, substantial progress has been made, both 
experimentally and theoretically, towards the understanding of the behavior of hard-sphere-like 
fluids ^ ^ as well as of dipolar mixtures under confinement.-"^ 

Confined systems comprised of soft particles have received substantially less attention. -SiM 
Even in the cases in which soft potentials have been considered, a simplifying assumption has 
been made, namely that the confined entities are point particles whose interaction potential does 
not depend on geometric confinement. However, for a wide class of common soft particles, ranging 
from polymer brushes to dendrimers and from microemulsions to microgels, the confining walls 
physically deform the particle, resulting into conformational changes that affect their interaction in 
substantial ways. No attempt has been undertaken to-date to take this important aspect into account 
and to then proceed towards an investigation of the structure of soft, deformable particles in the 
borderline between two and three dimensions. At the same time, soft systems possess important 
unique properties due to their inherent tunability. Thus, the effective interaction potential between 
two star polymers is known to change significantly in going from bulk (three dimensions) to planar 
confinement (two dimensions).— Therefore, this tunability presents both advantages and additional 
challenges. 



The central goal of the present work is to study structural properties of star polymers under 
confinement, where the confined conditions are created through adsorption of the star polymers on 
a planar wall, with the degree of confinement controlled via the star- wall adsorption strength. The 
behavior of star polymers in the vicinity of a substrate is important for a variety of scientific and 
technological applications. For example, formation of colloidal aggregates with desired ordered 
structure on a substrate can be achieved by controlling effective colloidal interactions via the use 
of star polymers which act as depletants. As such, adsorption of star polymers has been actively 
studied both experimentally ^^^^^ and theoretically.—"— 

The present study is motivated by the recent work that has demonstrated a dramatic amplifi- 
cation of tension along the backbone of branched macromolecules as a result of their adsorption 
on a substrate.—"— By analogy, one would expect that adsorption on a surface would result in an 
enhancement of repulsive interaction between the two star polymers. — In the work at hand, this 
conjecture is confirmed and quantified both via Molecular Dynamics (MD) simulations and scal- 
ing theoretical calculations. In particular, based on scaling arguments, we propose a simple and 
accurate functional form for describing mean force between two adsorbed stars, which is found to 
be in excellent agreement with the simulation results for all the combinations of parameter values 
studied. 

Computational model 

We consider a three-dimensional coarse-grained model— of a polymer star which consists of / 
linear arms with one end free and the other one tethered to a microscopic core (seed monomer) of 
size Re, which is of the order of the monomer size o. Each arm is composed of A^ particles of equal 
size and mass, connected by bonds. Thus, the total number of monomers in the star (excluding the 
core) is fN. Previous studies on star polymers— has shown that this model capture all the charac- 
teristics of the stars' conformations and interactions. The bonded interaction between subsequent 
beads is described by frequently used Kremer-Grest model— V'^°{r) = V^^'^^{r) + y^'^'^(r), with the 



so-called 'finitely-extensible nonlinear elastic' (FENE) potential: 

y--(r) = -0.5^r2ln[l-(r/ro)2]. (1) 

The non-bonded interactions between monomers are taken into account by means of the Weeks- 
Chandler-Anderson (WCA) interaction, i.e., the shifted and truncated repulsive branch of the 
Lennard- Jones potential, given by: 
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{O /rY^ - {O /rf + \/A\e{7}I^O -r). (2) 



In Eqs. (??) and (??), r denotes the distance between the centers of two monomers (beads), 
while the energy scale e and the length scale o are chosen as the units of energy and length, 
respectively. Accordingly, the remaining parameters are fixed at the values k = 30e/a^ and 
ro = 1.5. In Eq. (??) we have introduced the Heaviside step function 0(x) = or 1 for x < 
or X > 0. In consequence, the steric interactions in our model correspond to good solvent condi- 
tions. The substrate in the present investigation is considered simply as a structureless adsorbing 
plane, with a Lennard- Jones potential acting with strength Es in the perpendicular z-direction, 
V"(z) = 4£s[{o / zY^ — (a/z)^]. Since in the the majority of our simulations we considered the 
case of strong adsorption, £s/{kBT) = 5 with kg = I and T = O.lle/k^, the configurations of the 
simulated stars were two-dimensional. The chosen value e^ resulting into an adsorption energy of 
5k^T corresponds to values in the typical experimental range.— For comparison and complete- 
ness, we have also examined the case of stars confined between two parallel, repulsive walls (slit). 
In this case, stars form quasi two-dimensional configurations, since small lateral excursions of the 
order of the slit width H are allowed. The wall-monomer was modeled by WCA-potentials given 
by Eq. ?? acting only in z-direction and the slit width was set to H = 4<7. 

The equilibrium dynamics of the chains is obtained by solving the Langevin equation of motion 
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Figure 1 : Snapshots of thermalized star-polymers with / = 9 and N = 15 adsorbed on a solid 
substrate displayed for different separation distances: (a) d = 35a; (b) d = 3a. 



for the position r„ — [xn,yn,Zn] of each bead in the star, 



mr„ 



=+f: 



Crn + Rn{t),n=l,...JN 



(3) 



which describes the Brownian motion of a set of interacting monomers. Even though Molecular 
Dynamics would be equally well-suited for gathering equilibrium statistics, the approach based 
on the Langevin equation allows for the introduction of a fixed temperature via the coupling to a 
random force (see below) in a natural way. In Eq. (??) above, Y™^ and F*'^* are the deterministic 
forces exerted on monomer n by the remaining bonded and nonbonded monomers, respectively. 
The influence of the solvent is split into a slowly evolving viscous force — ^r^ and a rapidly fluc- 
tuating stochastic force. This random, Gaussian force R,, is related to friction coefficient C, by 
the fluctuation-dissipation theorem. The integration step employed was At = 0.002 time units. 



whereby time is measured in units of ^/mo^/e, m= \ denoting the mass of the beads. The ratio 
of the inertial forces over the friction forces in Eq. (??) is characterized by the Reynolds num- 
ber Re = \fm£.IC,<3 which in our setup is Re = 4. In the course of simulation the velocity- Verlet 
algorithm was employed to integrate the equations of motion (??). 

Starting configurations were generated as radially straight chains fixed to the immobile seed 
monomers. Stars were placed at a distance z = a above the surface and at large separation dis- 



tance between seeds d, in order to avoid monomer overlaps between neighboring molecules. As a 
next step star-polymers were equilibrated until they adopted their equilibrium configurations, from 
which sizes and monomer profiles have been determined. The required equilibration time depend 
on / and A^^, i.e., 5 x 10^ time steps for the smallest stars and up to 2 x 10^ for the largest. As a next 
step star-polymers were equilibrated until they adopted their equilibrium configurations. This step 
was controlled by monitoring the time evolution of the center-to-end radius R^ and the gyration 
radius Rg, as well as of the and radial density profiles of the stars. 

Well-equilibrated configurations from a large interstar separation were used as starting points 
for the next simulation, for which we move the second seed by distance AJ closer to the first one 
and repeated the equilibration step. In order to sample the interaction force for each equilibrated 
configuration (at a given distance d) a production run was performed during which time averaged 
force {F{t,d)) = F{d) exerted on the seed was computed. To assure that our measured forces 
were true equilibrium results, we have also employed the reverse path of first bringing the two 
seeds very close to one another and then increasing their mutual distance. The resulting force vs.- 
distance curves turned out to be independent of the procedure of moving the stars close or away 
from each other, excluding thus any hysteresis effects in our simulations. 

During the force measurement all particles (including the cores of the stars) were allowed to 
fluctuate in the z-direction. The fluctuations of seeded monomers in the z-direction do not con- 
tribute to the mean force F{d), since the time average of the z-component of the force, {F^it., d)), 
for seeded monomers vanishes. The run was continued until the uncertainty of the average forces 
for all three components dropped below O.OOSe/o. Our checks have also shown that time aver- 
aged net force acting on the seed is equivalent to the time averaged force felt by the core of a star 
polymer. Typical simulation snapshots of two interacting stars are presented in[T] 



Theory and comparison with simulations 

The theoretical approach we adopt for describing both single-star conformational properties and 
the effective interaction between two strongly adsorbed stars is based on the Daoud-Cotton blob 
picture^ of multiarm stars and its generalization to two interacting stars,— appropriately adopted 
to two spatial dimensions. Each of the / arms around the center of a star can be considered as being 
built of a succession of concentric blobs of radius ^ (s) , where s denotes the distance from the star 
center. The polymer segment confined within each blob is performing a two-dimensional, self- 
avoiding random walk (SARW), whereas the blobs are space-filling within the star corona. The 
latter assumption leads to the scaling ^ (s) = s/f; the former implies that the number of monomers, 
g{s), from a single chain at distance s from the center, scale as g{s) ~ ^ ^'^{s), where v = 3/4 is 
the SARW-critical exponent in two dimensions. Concomitantly, the monomer concentration c{s) 
scales with the distance from the star center as: 

C(s) ~ lil!ll^ ^ r2-l/v l/v-2 -I/V .4s 

where a denotes the segment size. Using the value of the scaling exponent v=3/4 in two dimensions 
under good solvent conditions, ^- one obtains: 

c{s) ^ f/' {^f\-\ (5) 

The size of the star R, standing for, e.g., its gyration radius Rg or its center-to-end radius Re, can 
be determined by the conservation law of the monomer number, fN = 2n Jq c{s)sds. The scaling 
of 7? follows: 

Rr^afl^N^I\ (6) 

featuring a stronger stretching effect of the arms with / in comparison to the three-dimensional 
case,— in which R ~ af^'^N^'^ . The above relations for c{s) and R imply that the scaled monomer 



concentration R^c{s)/{Nf) is a universal function p{x) of jc = r/R, with the properties: 



2k / p(x)xdx = 1 and p(x) ~ jc ^/^ for;c < 1. 
Jo 



(7) 



Our MD simulation results in |2] fully confirm this picture, since density profiles from three 
different stars collapse onto one another when properly scaled, showing in addition the power- 
law predicted by Eq. (??). Furthermore, the validity of Eq. (??) is demonstrated in|3l where the 
quantity f~^'^Rg is shown to indeed scale as N^'"^. Scaling- and blob-based considerations are 
thus fully capable of predicting quantitatively the properties of strongly adsorbed, individual star 
polymers. 
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Figure 2: MD simulation results for the normalized monomer density profiles for various star 
polymers as a function of dimensionless distance from the star center; the dashed line indicates the 
predicted power-law, c{s) ~ s^^'^ . 



The theoretical approach to the functional form of the effective interaction jSVeffl^) between 
two star polymers in two dimensions, whose anchoring points are separated by distance d, proceeds 
along two parallel lines, one pertaining to the case / = 1 (two polymer chains) and the second to the 
case / ^ 1 (multiarm stars). For / = 1 , we consider first the case of a single, self-avoiding chain of 
A'^ segments, which possesses a partition function Z\ [N) that scales as Z\ (N) ~ T^N'^^^ M Here, T 
is a model-dependent parameter and 7 = 43/32 is another universal exponent of the 2d-SARW.— 
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Figure 3 : Simulation results for the star radius R^ scaled by / ' as a function of the star arm length 
A'^ for several values of / and N\ the dashed line indicates the power law of Eq. (??). 



Two such stars, when kept at infinite separation, are thus characterized by a partition function 
T^iid — > oo) = Z-^{N) ~ x^^ N^^"^ \ When their end-monomers are kept at separation r r^ a, the 
two chains resemble a single SARW of IN steps, thus Z2{d ~ a) ~ x^^ {2Ny^^ . On the other hand, 
the segment length can be renormalized to larger values as long as it does not exceed, roughly, the 
equilibrium chain size R. The effective interaction is given as jSVeffl*^) = — ln[Z2(J)/Z2(<i -> °°)]. 
Using the above considerations in conjunction with the relation R ~ N^ , we obtain 



j8yeff(^) = -^ln(^) fox d<R, 



(8) 



for the case /= 1, where (7— l)/v = 11/24. 

In the opposite case, / ^ 1, we proceed by first counting the number of Daoud-Cotton blobs, 
M{f), of a single star. The condition c{s) ~ ^^'^{s) together with the geometrical picture of 
closely-packed blobs lead, after some straightforward algebra, to the result M{f) ~ f^\n{R/Rc), 
where Re is a (microscopic) scale denoting the size of the core of the star, on which the chains 
are grafted; this result should be compared with the corresponding scaling M{f) ~ f^'^\n{R/Rc) 
in three dimensions,=2i demonstrating, once more, the stronger stretching of the chains in two 
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dimensions. If two such multiarm stars are now brought to a separation of order Re from one 
another, they resemble a single star with 2/-arms, forming a number M{2f) ~ {2f)^\n{R/Rc) 
blobs. The effective interaction can be estimated by assigning the usual free energy cost of k^T 
per blob, i.e., pVeff{d) ~ M{2f) — 2M{f), leading to the expression: 

/3\/eff(rf) = -©(/) In (^) for d<R, (9) 

where the amplitude scales as ©(/) ~ /^ for / ^ 1. 

Since a single chain (/ = 1) can be seen as a special case of a star polymer and the short-range 
behavior of the effective potential is logarithmic with r in both cases, it is reasonable to request that 
0(/ = 1 ) = 1 1 /24 and introduce a general expression for this amplitude that smoothly interpolates 
between the two extreme cases. Field-theoretical arguments lead to the expression:— 

©(/) = -^, (10) 

which, evidently, satisfies the requirements on both limits, f —^ I and / ^ 1 . The ultrasoft logarith- 
mic repulsion for the star-star interaction potential has been obtained earlier in three dimensions on 
a similar basis, with the corresponding short-range repulsive force behaving as (5/^' ^/ 1 8) \n{R/d) .— 
The functional form of the short range star-star interaction remains the same but its amplitude in- 
creases due to confinement from 5/^/^/18 to (2 -|- 9/^)/24.— 

The effective force Fshort(^) for J < i? follows as minus the gradient of Veffl*^)- Taking into 
account that for d ^ Re the strong steric interaction between the cores dominates over the universal, 
logarithmic form of the latter arising from polymer self-avoidance, we are led to introduce a shift 
d ^ d — Re on the coordinate, as previously done also for three-dimensional star polymers,— 
yielding: 

kgT 24(J-i?,)- 

In the long-range regime, Flory theory yields a Gaussian form for the potential of mean force. 
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which resuks into the following scaling expression for the mean force at long {d > Rg) star-star 

^long(^)_ 2 + 9/2 d 



separations:— 



■exp 



{d/Ref 



(12) 



keT 24{R,-R,)R, 

Here, the prefactor has been determined using the physical requirement of continuity between the 
short- and long-range forces at d = Re. The corresponding results for the mean force between two 
adsorbed stars are shown in |4] together with the simulation data for several values of / and A'^; the 
agreement between theory and simulation is excellent for all values of parameters studied, without 
use of fit parameters. 

Setting X = d/Rg and taking the limit Re —^ (long chains), we put forward a simple and 
accurate analytical form for the effective potential between two star polymers in two dimensions 



as: 



2-uQf2 -lnx + e"V2, ifjc<l, 
i8Veff(x) = -^ I (13) 

e--'^V2, ifx>l. 

Eq. (??) extends the hitherto known three-dimensional result— to two dimensions, introducing 
an easy to handle functional form for the latter. The ultrasoft potential of a logarithmic form is a 
universal feature also for the interaction of dissimilar star polymers— or between star polymers and 
flat or curved walls— in three dimensions. The criterion of a 'weak mean field fluid',— expressed 
through the condition V = jQX^~^Veff{x)dx < oo, is fulfilled for the effective interaction of Eq. 
(??). 

Finally, in [5] we compare our simulation results for the two distinct cases already mentioned 
above: stars strongly adsorbed on a single wall and stars trapped between two parallel, repulsive 
walls forming a slit of width H = Ao in the z-direction. The data indicate a systematic deviation of 
the two forces, with growing inter-star separation distance d. As long as the two stars are far apart 
in the slit, the arms on their periphery experience very low monomer density in their vicinity, and 
can easily slip above or beneath one another. This is reflected by the smaller force of repulsion, felt 
by the stars, compared to the case of strong adsorption. In the latter case, the arms cannot intersect 
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Figure 4: MD simulation (symbols) and scaling theory (lines) results for the mean force between 
two adsorbed stars for several values of / and A'^, as indicated in the legends. 




Figure 5: Comparsion of MD simulation results for the mean force between two stars adsorbed 
onto surface (circles) or trapped in a slit (solid line). Here / = 9 and A^ = 50. Slit width is 4a. 
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as they stick strongly to the substrate, so they either have to bend when encountering an arm of the 
second star (meaning stronger repulsion), or crawl over it which is less probable because it costs 
energy of contacts with the substrate. On the other hand, the forces coincide at smaller distances 
d because the increased density around the star centers largely prevents easy threading of arms 
through each other, so that there is no way to avoid the strong repulsion. 

Conclusion 

We have demonstrated by a combination of molecular simulation and analytical theory that soft 
colloids, represented here by the common example of star polymers, experience dramatic confor- 
mational changes upon confinement in two dimensions, which profoundly affect their effective 
interaction, rendering it significantly more repulsive than in three dimensions. Previous studies on 
three-dimensional star polymers— indicated that the pair forces continue to provide a valid and 
accurate description of the interactions even for the case in which three stars approach each other 
at distances that represent a small fraction of their respective sizes. Thus, our results for effectively 
two-dimensional star polymers are expected to be quantitatively reliable for concentrations well 
above the overlap density of the stars. In addition our checks have shown that effective interaction 
between two adosorbed stars is the same as the effective force for stars trapped in a very narrow slit 
formed by two paralell and repulsive walls. Accordingly, we expect that much stronger ordering 
will take place in confinement and that the structural arrest lines will shift to lower concentrations 
with respect to their bulk values. Changing the amplitude of the attractive interaction to the wall 
allows us to tune the conformations and the interactions of the (partially) adsorbed stars, opening 
a pathway to tunable patterning of surfaces with ordered arrays of adsorbed stars. 

Possible experimental validation of our threoretical predictions would be to graft long chains 
on colloidal particles and use optical tweezers in order to measure effective force. The alternative 
approach would be to measure the correlations between the adsorbed stars in the same way this 
has been done for other 2d-colloidal systems adsorbed on surfaces or interfaces (see e.g.—). The 
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resulting correlations should then be describable by the effective interactions derived in this paper. 
Our work opens up a possibility to consider also the case of partially absorbed stars with fa 
chains lying on the plain and / — fa chains assuming 3d-configurations, for which the force can be 
obtained as the sum of the 2d- and 3d-expressions, each weighted with the appropriate amplitude. 
The determination of the interaction between partially adsorbed stars, bridging the gap between 
two and three dimensions, will be the subject for future work. 
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